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Abstract. 

Thanks to the high multiplex and efficiency of Giraffe at the VLT we have been able for the 
first time to observe the Li I doublet in the Main Sequence stars of a Globular Cluster. At the 
same time we observed Li in a sample of Sub-Giant stars of the same B-V colour. 

Our final sample is composed of 84 SG stars and 79 MS stars. In spite of the fact that SG and 
MS span the same temperature range we find that the equivalent widths of the Li I doublet in SG 
stars are systematically larger than those in MS stars, suggesting a higher Li content among SG 
stars. This is confirmed by our quantitative analysis carried out making use of ID hydrostatic 
plane-parallel models and 3D hydrodynamical simulations of the stellar atmospheres. 

We derived the effective temperatures of stars in our the sample from Ha fitting. Theoretical 
profiles were computed using 3D hydrodynamical simulations and ID ATLAS models. Therefore, 
we are able to determined ID and 3D-based effective temperatures. We then infer Li abundances 
taking into account non-local thermodynamical equilibrium effects when using both ID and 3D 
models. 

We find that SG stars have a mean Li abundance higher by 0.1 dex than MS stars. This result 
is obtained using both ID and 3D models. We also detect a positive slope of Li abundance 
with effective temperature, the higher the temperature the higher the Li abundance, both for 
SG and MS stars, although the slope is slightly steeper for MS stars. These results provide an 
unambiguous evidence that the Li abundance changes with evolutionary status. 

The physical mechanisms responsible for this behaviour are not yet clear, and none of the 
existing models seems to describe accurately these observations. Based on these conclusions, we 
believe that the cosmological lithium problem still remains an open question. 
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f Present address: Dpto. de Astrofisica y Ciencias de la Atmosfera, Facultad de Fisica, Uni- 
versidad Complutense de Madrid, E-28040 Madrid, Spain. Email: jonay@astrax.fis.ucm.es 



2 



Gonzalez Hernandez et al. 



1. Introduction 

The determination of the baryonic density from the fluctuations of the cosmic mi- 
crowave background (CMB) by the WMAP satellite (Spergel et al. 2007, Cyburt et al. 
2008) implies a primordial Li abundance which is log(Li/H) + 12 = 2.72 ± 0.06, at least 
0.3-0.5 dex higher than the Li abundance determined in metal-poor stars of the Galactic 
halo (Spite & Spite 1982). 

Many different models of Li depletion have been proposed to explain discrepancy: 



(a) Piau et al.(2006) proposed that the first generation of stars. Population III stars. 



could have processed some fraction of the halo gas, lowering the lithium abundance; 
(b) other authors suggest that the primordial Li abundance has been uniformly de- 
pleted in the atmospheres of metal-poor dwarfs by some physical mechanism (e.g. tur- 
bulent diffusion as in Richard et al.(2005)[ Korn et al.(2006) gravitational waves as in 
Charbonnel fc Talon(2005)[ etc.); and (c) finally, it has been also suggested that the stan- 
dard Big Bang nucleosynthesis (SBBN) calculations should be revised, possibly with the 
introduction of new physics as in e.g. Jedamzik(2004) Jedamzik(2006) Jittoh et al.(2008) 



Hisano et al.(2009) 



Here we present the determination of Li abundances of subgiant (SG) and main- 
sequence (MS) stars of the cluster NGC 6397. This work provides the first observations 
of the Li doublet in MS stars of a Globular Cluster. 



2. Observations 

We performed spectroscopic observations of the globular cluster NGC 6397 with the 
multi-object spectrograph FLAMES-GIRAFFE at the VLT on 2007 April, May, June 
and July, covering the spectral range AA6400-6800 A at resolving power X/5X ~ 17,000. 

We selected subgiant and dwarf stars in the colour range B — V = 0.60 ± 0.03, which 
ensures that both set of stars fall in a similar and narrow effective temperature range 
(see Fig. 3 online in Gonzalez Hernandez et al. 2009a). 

The spectra were reduced with the ESO pipehne and later on treated within MIDAS. 
We correct the spectra for sky lines and, barycentric and radial velocity. We typically 
combine 17 spectra of dwarfs and 4 spectra of subgiants to achieve a similar S/N ratio 
between 80 and 130 in both sets of stars. The mean radial velocity of the cluster stars is 
Vr = 18.5 km s^^. 



3. Stellar parameters 

We derived the effective temperature by fitting the observed Ha line profiles with 
synthetic profiles, using 3D hydrodynamical model atmospheres computed with the 



CO^BOLD code. The details of this code are provided in^Freytag et al.(2002) and Wedemeyer et al.(2004) 



The ability of 3D models to reproduce Balmer line profiles has been shown in Behara et al.(2009) 



In that work, the Ha profiles of the Sun as well as the metal-poor stars HD 84937, HD 
74000 and HD 140283 are studied. Ludwig et al.(2009) also quantified, from a purely 
theoretical point of view, the difference between the effective temperatures derived from 
Ha fitting using ID and 3D models. 

We also derived the effective temperatures of MS and SG stars using ID ATLAS 9 
model atmospheres (see Kurucz 2005) and the same fitting procedure. In Fig. [1] we dis- 
play the histograms of the effective temperatures derived for MS and SG stars using both 
ID and 3D models. In the ID case we got similar effective temperatures for both sets 
of stars. However, using 3D hydrodynamical models we obtained hotter temperatures by 
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Figure 1. Histograms of ID (left panel) and 3D (right panel) effective temperatures derived by 
fitting the observed Hq profiles with theoretical profiles computed using f D (left panel) and 3D 
(right panel) model atmospheres, in bins of fOO K, in MS (solid line) and SG (dashed-dotted 
line) stars in the Globular cluster NGC 6397. 

approximately 250 K in MS stars and 150 K for SG stars. In the 3D case the histogram 
of SG stars is slightly shifted with respect to the histogram of MS stars to cooler tem- 
peratures, as expected from the difference in surface gravity between MS and SG stars 
and the sensitivity of the B — V colour to the surface gravity. 

Fixed values for the surface gravity were adopted for both subgiant and dwarf stars 
in the sample, according to the values that best match the position of the stars on a 12 
Gyr isochrone (Straniero et al. 1997). The adopted values were log((7/cm s^) ~ 4.40 and 
3.85 for MS and SG stars, respectively. 



4. Li abundances 

We measure the equivalent width (EW) of the Li i 6708 A line in SG and MS stars by fit- 



ting synthetic spectra of known EW to the observed Li profiles. Gonzalez Hernandez et al. (2009a) 



showed the histograms of the EWs measured in SG and MS stars of this cluster (see their 
Fig. 1). In that figure it is clearly seen that the EWs of SG stars are larger than those 
of MS stars. They also estimate the weighted mean EW of the SG stars, being ~ 1.1 pm 
larger than the weighted mean EW of MS stars. Although the colour i? — is sensitive to 
surface gravity, a priori, this result was not expected, and clearly suggests that subgiants 
in this cluster have actually higher Li abundances than dwarfs. 

We derived Li abundances using 3D model atmospheres. The line formation of Li 
was treated in non-local thermodynamical equilibrium (NLTE) using the same code and 
model atom used in iCayrel et al.(2007) The analysis was also done using ID model 
atmospheres providing essentially the same picture, even when T^ff in ID are lower (see 
also Fig. 6 online in Gonzalez Hernandez et al. 2009a). In the ID case we used the 



Carlsson et al.(1994) NLTE corrections. 



5. Discussion and Conclusions 

In Fig. [2] we display the histograms of the derived ID-NLTE and 3D- NLTE Li abun- 
dances of dwarf and subgiant stars of the globular cluster NGC 6397. In both ID and 3D 
cases, the SG stars have on average larger amounts of Li content in their atmospheres 
than MS stars. The difference in the mean Li abundance of dwarfs and subgiants is 
~ 0.14 dex in the ID case and ^ 0.07 dex in the 3D case. This difference between the ID 



4 



Gonzalez Hernandez et al. 




Figure 2. Histograms of ID (left panel) and 3D (right panel) non-LTE Li abundances, in bins 
of 0.07 dex, in MS (solid line) and SG (dashed-dotted line) stars in the Globular cluster NGC 
6397 



and the 3D cases is due to the effective temperatures derived using ID and 3D models, 
and is not related to the NLTE corrections which were independently computed in the 
ID and 3D cases. In fact, due to the cooler temperatures derived in ID with respect to 
the 3D case, the mean ID Li abundance lower by ^ 0.13 dex for the MS stars ^ 0.06 dex 
for the SG stars than the mean 3D abundance. 



Lind et al.(2009) also find different mean Li abundances, using ID models, in MSs 



and SGs, but only by 0.03 dex although still significant at la. However, their result is 
partially affected by the very narrow range of Toff for MSs deduced by iLind et al.(2009) 



(~ 80 K) compared to the wide range {'-^ 450 K) for the SGs (see Fig. 7 online in Gonzalez 
Hernandez et al. 2009a). 



Gonzalez Hernandez et al. (2009a) showed, in their Fig. 2, the 3D-NLTE abundances of 



SG and MS stars versus 3D effective temperatures (see Fig. 1 in Gonzalez Hernandez et al. (2009b) 



for a similar picture but with Toff and Li abundances computed using ID models). The 
points in that figure display a decreasing trend of Li abundance with decreasing tem- 
perature. This lithium abundance pattern is different from what is found among field 
stars (see e.g. Melendez & Ramirez 2004, Bonifacio et al. 2007, Gonzalez Hernandez et 
al. 2008). 

Our results imply that the Li surface abundance depends on the evolutionary status 
of the star. In Fig. 2 of Gonzalez Hernandez et al. (2009a) the Li isochrones for different 
turbulent diffusion models (Richard et al. 2005). These models were shifted up by 0.14 
dex in Li abundance to make the initial abundance of the models, log(Li/H) = 2.58, 
coincide with the primordial Li abundance predicted from fluctuations of the microwave 
background measured by the WMAP satellite (Cyburt et al. 2008). The models assuming 
pure atomic diffusion, and, among those including low levels of turbulent mixing, in 



particular the model T6.0, preferred model in Korn et al.(2006) and Lind et al.(2009) 



are ruled out by our observations, since they predict the Li content in MS stars to be 
higher than in SG stars. The only model that predicts a Li pattern which is qualitatively 
similar to that observed, is the T6.25 model. This model shows a decreasing trend of 
Li abundances with decreasing temperatures and also predicts higher abundances for 
SG stars than for MS stars although at temperatures lower than 6000 K, which is not 
consistent with the observed trend even in the ID case (see Fig. 1 in Gonzalez Hernandez 
et al. 2009b). 

Models including atomic diffusion and tachocline mixing (Piau 2008) do not seem to 
reproduce our observations either, since they provide a constant Li abundance up to 
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5500 K. More sophisticated models are required, for instance, those models that besides 
including diffusion and rotation also take into account the effect of internal gravity waves 

(Talon & Charbonnel 2004) , seem to predict accurately the Li abundance pattern in solar- 
type stars, at solar metallicity (Charbonnel & Talon 2005), but models at low metallicity 
are still needed. 

The cosmological lithium discrepancy still needs to be solved. Given that none of the 
existing models of Li evolution in stellar atmospheres matches our observations in the 
globular cluster NGC 6397, we hope our results will prompt further new theoretical 
investigations. 
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